function y = phi1phi1(x)
    y = phi1(x) .* phi1(x);
end